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Abstract 

Integra! transform techniques such as the Laplace transform, provide simple and direct 
methods for solving viscoelastic problems formula tec within a context of linear materia; response 
and using ..rear measures for deformation Application of me transform operator reduces the 
governing linear integ* c~ d:f f erential ecuat<cns to a set of algebraic relat one between the 
transforms o' the unknown functions the v.scoclastic operators, and the initial and boundary 
conditions ’-aversion eitr.e- erectly or through the use 01 the appropriate convolution theorem, 
provides the time domain response once the unknown functions have been expressed in terms of 
sums products o- r £t>os o' known transforms When exact inversion >s net possible approximate 
techniques such as suggested by £ draper y may provide accurate results 

r ne over all problem Becomes substantially more comp-ex when noniirear effects must be 
included We consider here situations where a linear material constitutive saw can still be 
productively employed tut where the magnitude of trie resulting time dependent deformations 
warra-ts the use of a nonlinear Kinematic analysis The governing equations w.li be nonlinear 
mtegro- differential equations for this class of problems Thus traditional as well as approximate 
techniques such a; cited above cannot be employed since tne transform of a nonlinear function is 
not expncfiy expressible 

Rogers and lee : considered such a problem in an investigation of the finite deflection of a 
viscoelastic cantilever beam Employing an analogy to an associated elastic problem they derived a 
solution to the viscoelastic problem n a fo r m involving a time convolution of a nonlinear space and 
t me Dependent integral funct on. Numerical evaluation was accomplished using Picard's method of 
successive substitutions. Newton- Cotes quadratures \ve r e employed to approximate the spatially 
dependent integral relationship; a mean value based finite difference formula was used for the time 
convolution 

Solution procedures of this type are generally well so ted for computer implementation, 
however they can become computationally inefficient when the response must be determined over 
r i extended tune period Each increment in time requires a ree valuation of the convolution integrals, 
"’vs •vie entire deformation history must be retained in memory during the calculations. Since each 
completed se* of computations adds another set of results to this history this results in an ever 
inert a sing memory requirement, in addition the total number of computations which must be 
per : orr.es during the succeeding iteration is aiso increased. 

n this regard an approximation technique proposed by Scnapery x can provide an attractive 
alternative Commonly referred to as the ' quasi-elastic" approximat'd • it has most recently been 
employed by Vinogradov**’ and Vinogradov and V 'gewesra' ; .n studies of the behavior of 
eccentrically loaded viscoelastic cantilever beams 

T he method is based on the observation that the solution procedure developed by Rogers and 
.ex- may bo interpreted as a sequence of short time interval ' quasi-eiactic' solutions This suggests 
Thai approximate solutions may be generated by replacing the original viscoelastic problem by an 
equivalent' time dependent elastic one !n this replacement problem the elastic properties are 
equated to the instantaneous values of the relaxation moduli or creep compliances of the 
viscoelastic material 
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The inherent numerical advantage provided by tms technique is that it can eliminate the 
potentially inefficient convolution integral calculations Thus the speed and efficiency at which the 
time dependent response is calculated vvil' become independent of elapsed time The obvious 
potentia 1 disadvantage of the technique is that since it is an approximation significant differences 
may exist between the actual response of the viscoelastic body and those predicted 
quast-elasncally. In addition, the quasi-elastic method does not provide a direct method for 
assessing whether or not the errors which may be incurred are conservative. 

In this paper we demonstrate that the quasi-elastic approximation technique can be modified 
to prov:ae both upper and lower bound predictions for the class of viscoelastic problems under 
con«'deraticn lo accomplish this the solution of tne actual viscoelastic problem is first formicated 
m terms of a Volterra type integral equation of the second kina Ar. upper or lower bound solution 
is then established by replacing the desired unknown function which appears under the integral sign 
by an appropriate approximating funct'on. The approximating function selected is one which 
inherently bounds the aesired unknown function and is also independent of the integration variable 
The kerne! of the integral operator can then be '.ntegrated formally casting the approximate problem 
into the typical quasi-elastic format Provided the Kernel of the original .ntegral operator is positive 
semi-definite this quasi-elastic solution must bound the actual solution in tne same manner that the 
approximating function bounded the desired unknown function 

Illustrated on the following page are comparisons o f results obtained using the bounding 
technique to the exact solution for the transversely loaded viscoelastic cantilever beam problem 
considered by Rogers ana Lee Additional examples of the application of the technique, including 
procedures to be used when nonlinear boundary terms appear in the original .ntegral equation, are 
provided in the paper. 
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Fig 1 Viscoelastic response at midspan of Fig 2 Viscoelast.c response at end point 

transversely loaded cantilever c f transversely loaded cantilever 


